Local impurity effects in superconducting graphene 
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We study the effect of impurities in superconducting graphene and discuss their influence on 
the local electronic properties. In particular, we consider the case of magnetic and non-magnetic 
impurities being either strongly localized or acting as a potential averaged over one unit cell. The 
spin dependent local density of states is calculated and possibilities for visualizing impurities by 
means of scanning tunneling experiments is pointed out. A possibility of identifying magnetic 
scatters even by non spin-polarized scanning tunneling spectroscopy is explained. 



Graphene - one monolayer of carbon atoms arranged 
in a honeycomb lattice - can be considered as one of 
the most promising materials for post-silicon electronics 
[1] . Impurity states often are the basis of modern elec- 
tronic devices. Because of its peculiar electronic struc- 
ture graphene provides a model system for studying rel- 
ativistic quantum physics in condensed matter [2, 3]. A 
careful investigation of the impurity states in graphene 
also allows one to address fundamental questions of rela- 
tivistic scattering [4-8]. Thus, far the central component 
to these studies is graphene's linearly vanishing density 
of states (DOS) near the Dirac point: It allows for virtual 
bound states (VBS) in this region, i.e. resonances in the 
DOS, which may be arbitrarily sharp in the vicinity of 
the Dirac point. The broad similarities of these impurity 
states with states in d-wave superconductors have been 
used to discuss common aspects of impurity scattering in 
superconductors and in graphene [5, 6, 8-10]. 

Recently, it has been demonstrated that it is possible 
to induce superconductivity in graphene via a proximity 
effect [11]. Prior analysis of the ballistic superconduct- 
ing transport [12] revealed an interesting suppression of 
the critical current near the Dirac point. It is there- 
fore simply a matter of time before the defects and im- 
purity induced states in superconducting graphene will 
be addressed locally. The case of superconductivity in 
graphene where opposite valleys arc nontrivially involved 
is also an explicit example of supcrcondcutivity in val- 
lyetronics [13]. The observation of the proximity effect 
in graphene raised fundamentally new questions about 
impurity effects in this material in presence of supercon- 
ducting correlations: i) Is there a possibility of intragap 
bound states and ii) what is the impact of gap opening 
on the Friedel oscillations in the continuum? 

In this letter we show that magnetic impurities do pro- 
duce impurity bound states inside the superconducting 
gap. These bound states always coexist with the for- 
merly studied VBS in the continuum. Thus the predicted 
impurity states are similar to the magnetic impurity in- 
duced states, so called Yu Lu-Shiba-Rusinov states, in 
s-wave superconductors [9, 14]. Due to its 2 dimension- 
ality graphene is well suited for Scanning Tunneling Mi- 



croscopy (STM) investigations and first experiments on 
normal state graphene already indicate the importance 
of impurity effects in this context [15, 16]. Therefore, we 
elucidate the real space shape of these impurity states, 
which will be directly observable in STM experiments. 

While graphene intrinsically is not superconducting, 
the Ti/Al bilayer contacts placed on the graphene sheet 
induce a measured supercurrent [11]. No spectral gap in 
the samples has been measured to date. We argue that 
the residual electron-electron interaction in the graphene 
will produce a gap in the spectrum. This gap will be pro- 
portional to the interaction strength and it remains to be 
seen how large this gap can be in the graphene. Electron 
spectroscopy such as STM and/or planar tunneling into 
graphene in proximity to superconducting leads would 
be able to reveal the spectroscopic gap. We will treat 
superconducting gap A below as a phenomenological pa- 
rameter that needs to be determined separately. 

Low energy electronic excitations in graphene can be 
described by two species of Dirac fermions located around 
two nodal points K in the Brillouin zone with the speed 
of light being replaced by the Fermi velocity Vf and the 
corresponding Hamiltonian H K ± = Vfh~(ki<Ji ^feca). °*> 
i = 1,2,3, are Pauli matrices acting on the sublattice 
degrees of freedom and o~a is the identity matrix. To 
understand impurities in superconducting graphene, we 
use the Nambu formalism including both valleys: 

H = -ihv f J d 2 x&{x)(di<Ji®T - d 2 o- 2 ®T 3 )®A 1 &(x) 

with $(a;)t = (^\ K+ (x),'fl K _(x),^^ K -{x),^^ K+ (x)) 
and ty^ K ±(x) being field operators of electrons with a 
spin fj, and belonging to a valley IT . T{ and A, with 
% = 1,2,3 are Pauli matrices but acting on the valley 
and Nambu space, respectively. To and Ao are the cor- 
responding identity matrices. In contact with a super- 
conductor, the proximity effect imposes a finite pairing 
potential A<73 Cg> To (8 Ai to the graphene sheet and re- 
sults in electron dynamics being described by the Dirac- 
Bogoliubov-de Gennes (DBdG) Hamiltonian [17]: 

H = -ihvf(dicri<&T -d20-2®T 3 )®Ao+Ao- 3 <g>To<g>A 1 . (1) 
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Figure 1: (Color online) Among the various local impurities 
we discuss two limiting cases. The scalar impurity (left), V s , 
corresponds to a uniform potential averaged over one unit cell, 
whereas the on-site impurity (right), V , acts on one sublattice 
only. 

To elucidate the effect of different impurities, we discuss 
both, a homogeneous potential acting within one unit 
cell V s (referred to as scalar impurity, Fig. l.a) as well 
as a strongly localized impurity V (referred to as on-site 
impurity, Fig. l.b) acting only at sublattice A and giving 
rise to intervalley scattering. 

Starting from impurity operators in the tight-binding 
form of e.g. Ref. [8] and using the conventions of Eqn. 



(1) we obtain the following explicit expressions for the 
impurity potentials in the adopted matrix notation: V s = 
V a ® to <8> A 3 + Viao <g> To ® A and V a = V Q (a 3 + a ) (g> 
(t +Ti) <g> A 3 + Vi(a 3 + oo) ® (r + n) <g> A [18] . In both 
cases Vo and V\ describe the electrostatic and magnetic 
contribution to the impurity potential, respectively. 

The effects of these impurities on the local electronic 
properties of the superconducting graphene sheets is con- 
tained in the local density of states (LDOS), which we 
calculate using the T-matrix approach [9]: In operator 
form, the full Green's function G(ui) in presence of the 
impurity is obtained from the unperturbed Green's func- 
tion G°(uj) via G(u) = G°(u) + G°(lu)T(lo)G°(uj) with 
T(uj) = V s{o) (1-G°(lu)V s{o) )- 1 . 

Dealing with local impurities, it is convenient to 
adopt the position space representation. Therefore, 
the free x-dependent Green's function G°(x,lu) in po- 
lar coordinates, x — x(r,<j>), is obtained from its mo- 
mentum space counterpart G°(p,uj) = (u> — H)^ 1 = 

(^tTQ^TQ+Df [pigl®To-p2<T2®T 3 ])®Ao + Ag3®To®Al , Found' 
LJ 2 — VfP 2 — A 2 ■> 

transformation 



G°( Xl L0) = 



g (r, w)(cjcro ( 8iTo(g)A -(-Acr3(g)ro(8)Ai)+5i(r, to) ([cos </>cri(g>T +sin (/xt 2 ®t 3 ]®A ) (2) 

I 



with g (r,uj) — v\ / Q Pc dppJ (pr)(W 2 (uj 2 — A 2 — vf p 2 )) 1 
and gi(r,uj) — ivf dpp 2 J\(pr)(W 2 (uj 2 — A 2 - 
vfp 2 ))~ l , where we expressed the Brillouin zone volume 
tt B = 2itW 2 /v 2 in terms of the bandwidth W. The 
Green's function at x — determines the LDOS of the 
free system and it occurs in the T-matrix: G°(0, u+iS) — 
M(aj){ua Q ® t q ® A + Atr 3 ® r ® Ai). Here is M(u>) = 
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A —it) 
W 2 +A 2 -u 2 



A 2 



and 
W 2 and 



M'(w) + iM"(w) with M'(w) = 

M"(w) = -^^^ for A 2 < lj 2 < 
M"(uj) = else. One sees, that the corresponding LDOS 
vanishes within the superconducting gap (ui 2 < A 2 ) and 
is given by Nq(lj) = outside the gap. 

In general, impurity resonances occur when the T ma- 
trix becomes (almost) singular, i.e. det(l — G°(0, oj)V) = 
0. For the scalar impurity this secular equation yields 
1 - 2M(lo)ujVi + M 2 (uj){lo 2 - A 2 ){V? - Vq 2 ) = with 
solutions that can be understood analytically in the fol- 
lowing limiting cases: Firstly consider a solely magnetic 
impurity, i.e. Vb = 0, with V\ > 0. In the Born limit the 
solutions loq = — A ± 6u> with 



5u> 



W 
2A' 



S -W 2 /(AV!) 



(3) 



the continuum approaching the gap edge exponentially 
with decreasing V\ . In the opposite limit of unitary scat- 
tering luq — ±A — 5u> with 



Suj = — 



2W Z 



(4) 



give rise to intragap bound and virtual bound states in 



fulfills the secular equation, where the upper (lower) 
sign corresponds to a intragap bound (continuum virtual 
bound) state. 

The numerical solutions to energies of the intragap 
bound states are shown in Fig. 2. It recovers the lim- 
iting cases obtained analytically and demonstrates also 
the effect of an electrostatic contribution Vb to the im- 
purity potential: In the Born limit, the exponential de- 
pendence of Sui on the magnetic potential strength V\ 
is dominant and suppresses any significant influence of 
Vq on the impurity state energy. In the V% — ► oo limit, 

Vb leads to a renormalization of the effective magnetic 

v 2 

potential strength Vi — > Vi(l — As Fig. 2 shows, 

the effect of an additional electrostatic potential becomes 
most pronounced in the intermediate region. There, the 
electrostatic contribution reduces the effective magnetic 
potential strength most significantly. 
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Figure 2: (Color online) Energy of the impurity resonance 
for the scalar impurity as a function of the magnetic impurity 
potential V\ f° r different electrostatic potentials Vq. The gap- 
parameter is A = W/1Q. The lower right inset shows the 
splitting of the impurity state due to intervalley scattering. 
We model intervalley scattering as V — Vi(ao ® to ® Ao + 
aao ® ri ® Ao) with the strength of the intervalley scattering 
parametrized by a and Vi = 5W. One split state is shifted to 
the gap edge, the other state remains an intragap state. 

Having understood the energy of intra gap bound 
states due to scalar impurities, we address now the 
strongly localized on-site impurities and elucidate the ef- 



fect of inter valley scattering. Due to valley degeneracy, 
the scalar impurity gives rise to doubly degenerate in- 
tra gap bound states. This degeneracy is lifted by in- 
tervalley scattering, see Fig. (2). The secular equation 
corresponding to the on-site impurity, 1 — 8M(u;)u;Vi + 
WM 2 (lu){u 2 - A 2 )(V? - V { 2 ) = 0, reduces to that of 
an on-site impurity with the replacement Vo,i — > 4Vo,i- 
Besides the lifting of the valley degeneracy, additional 
intervalley scattering results also in a renormalization of 
the effective impurity strength. 

With the real-space Green's function G(x, x' , u>) = 
G°(x - x', oj) + G°(x, lu)T{uj)G°{-x\ oj) one obtains the 
local density of states N(x,oj) = N q (uj) + SN(x,uj) = 
— ^ImG(x, x, oj) in presence of an impurity. This LDOS 
is a matrix corresponding to the matrix structure of the 
Green's function. It accounts for the contributions from 
the different sublattices, valleys and the Nambu space. 
According to the convention in Eqn. (1), the spin-up ex- 
citations are hole excitations yielding for each spin com- 
ponent the LDOS A^j (a:, uj) = Tr M^-N(x, ±uj), where 
the trace involves either the spin-down or up part of the 
Nambu space. In the case of the scalar impurity, this 
yields explicitly the following corrections to the unper- 
turbed LDOS in the continuum 



JAT iT (r,±c) = --Im- mu)uVi + M2(a;)(w2 _ A2)(v? _ vi) (5) 
with a iT = (u? - A 2 )[±Vu + M{uj)lu(V(? - Vf )] + {uj 2 + A 2 ^ and b SA = (±V + V{) + M{Lo)uj{V n 2 - V?). 



By replacing M(ui) — > 4M(w) in these formula, one ob- 
tains the case of the strongly localized on-site impu- 
rity. In STM experiments, graphene's lattice structure 
will give rise to a triangular modulation of the impurity 
states. This is neglected here, as similar effects in normal 
state graphcne have been discussed in Ref. [8]. 

Due to Eqn. (5), the asymptotic decay of the in- 
tra gap bound states at large distances from the im- 
purity is governed by go(r, ojq) and <? 2 (r, u)q). Neglect- 
ing high-energy cut-off related oscillations at this length 
scale, one may extend the momentum space integrals 
in Eqn. (2) to infinity. This yields modified Bessel 
functions, i.e. go(r,LUo) = — ykz~K(0, r-^/A 2 — u^/vf) and 

ffi(r,wo) = - V ^""° K(l,r^A 2 -^ 2 /v f ). Therefore 
the wavefunctions of the impurity states decay as 

6Nn(r,±u ) oc r -i e -W A3, -"S/»f. (6) 

As Fig. 3 (left) shows, impurity states in the gap give 
rise to prominent features in future STM experiments: 
The density of the impurity state at the impurity site at 



r = as well as the maximum of the density are strongly 
sensitive to the particular type of impurity. In general, 
impurity states with energies in the middle of the gap 
(Vi = 514^ in Fig. 3 (left)) give rise to the sharpest max- 
ima in the r-dependent LDOS. The ratio of the maximum 
density to the density at the impurity site increases with 
the potential strength V\. 

The ringstructure corresponding to these impurity 
states in STM images may give a powerful experimental 
tool for identifying particular impurities present in super- 
conducting graphene. This is in contrast to the normal 
state graphene, where weak impurities do not give rise to 
resonances near the Dirac point and will therefore hardly 
be apparent in scanning tunneling spectroscopy (STS)[8]. 

In the continuum, Eqn. (5) encodes the real space 
shape of VBS and Friedel oscillations around the im- 
purities. As Fig. 3 (right) shows exemplarily for a 
scalar impurity, the wavelength A of these oscillations 
is in any case determined by the energy uj and the gap 
A: A = nVf/y/uP — A 2 . Besides these oscillations giv- 
ing rise to standing wave patterns in future STM exper- 
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Figure 3: (Color online) Left panel: Density N of the intra- 
gap bound states as a function of the distance r from the 
impurity for purely magnetic scalar impurities and different 
potentials Vi . The impurity strength is given in units of the 
bandwidth W . Right panel: Friedel oscillations in the local 
density of states (LDOS) around a scalar impurity at r — 
with Vo = and Vi = 3W. The different curves correspond 
to the energies ui = 0.8, 0.5 and 0.25W. In both panels, the 
gap-parameter is A = IF/10. 
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the LDOS and LSDOS are not symmetric under particle 
hole transformation. Therefore, the degree of symmetry 
of the local spectra allows to estimate, whether the im- 
purity potential is magnetic or not - even in a non-spin 
polarized scanning tunneling spectroscopy experiment. 

In conclusion, we argued that magnetic scattering 
will produce impurity induced bound and virtual bound 
states in superconducting graphene. These impurity 
states are similar to the Yu Lu-Shiba-Rusinov states in 
s-wave superconductors [9] and exhibit an intricate real 
space and particle-hole dependent structure. We dis- 
cussed the energy dependence of these states as a func- 
tion of the potential parameters and pointed out char- 
acteristic oscillation pattern as well as decay properties 
in real space. This spectroscopic and topographic infor- 
mation can be obtained by STM [14]. We showed that 
each impurity generates a specific signature in the real 
space LDOS and provided a guideline for identifying dif- 
ferent impurities in future experiments. Since the Cooper 
pairs have zero momentum we find superconducting state 
in graphene to be a nontrivial example of valleytronics 
where valley quantum numbers are important [13]. 
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Figure 4: (Color online) The local density of states (LDOS, 
upper panel) and the local spin density of states (LSDOS, 
lower panel), 5N^(r,co) — 5Ni(r,co), at r = and r = 1 is 
shown for scalar of impurities with different potentials: a 
purely magnetic impurity with Vo = and Vi = 3W (left) 
as well as an impurity with Vo = 2W and Vi = 1W. 

iments, certain resonances due to VBS will be an even 
more prominent as well as impurity specific feature in 
these experiments: The LDOS in Fig. 3 (right) exhibits 
a characteristic peak at r ~ 1 and uo « 0.25VF. 

So far, we have discussed the r-dependent LDOS for 
different impurities and at different energies correspond- 
ing to STM images at fixed bias. The impurities will 
manifest themselves also in the energy dependence of the 
LDOS at fixed position, which is accessible by STS. In 
Fig. 4, the LDOS near a purely magnetic scalar impu- 
rity with Vq — and V\ = 3W is compared to an impu- 
rity contributing an electrostatic potential Vq = 2W and 
V\ = 1W. The purely magnetic impurity, Fig. 4 (left), 
does not break particle hole symmetry and yields there- 
fore a fully symmetric LDOS and a fully antisymmetric 
local spin density of states (LSDOS). This is in contrast 
to the more general second impurity, Fig. 4 (right), where 
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